Journal of The Royal Society Interface

Vaccine effectiveness (VE) estimates can vary widely between years and populations, even for the same vaccine. Estimated VE is known to be sensitive to susceptible depletion and differences in pre-vaccination infection risk between vaccinated and unvaccinated populations. However, how variation in pre-vaccination risk within and between the two groups affects VE estimates over time remains unclear. This uncertainty is especially important given negative VE estimates. We investigated the difference between estimated VE and true vaccine protection considering continuous distributions of pre-vaccination infection risk under three scenarios. When the vaccinated and unvaccinated populations differ in their mean risk, estimated VE can be higher or lower than true vaccine protection. Similar patterns arise when both populations share identical means but different risk distributions. Finally, if infection-derived immunity lasts longer than vaccine protection, annual VE estimates can vary by tens of percentage points between years despite constant true vaccine protection. These theoretical results underscore that VE studies estimate contrasting risk between vaccinated and unvaccinated individuals in a particular time and place, and VE estimates can vary counterintuitively between years and populations even with constant vaccine-induced protection. Explaining variability in estimated VE thus requires a more complete understanding of populations' distributions of infection risk.

The erythrocyte sedimentation rate (ESR) is one of the most common and widely used laboratory diagnostic parameters in connection with inflammatory reactions and it is probable that every reader has already experienced a determination of their ESR. A rapid ESR is a non-specific parameter that provides information about the inflammatory process. Although the origins of this methodology date back to antiquity, the description of the process as the collapse of a percolating gel formed from erythrocytes has only recently been achieved. It was not yet known whether slow ESR has any medically relevant significance. Here we show a variety of clinical pictures that exhibit a systematically slow ESR (e.g., sickle cell disease, neuroacanthocytosis syndromes, chronic mountain sickness). Using a combination of measured data and physical modelling, we show how the accuracy and significance of ESR data can be increased. With this improved ESR (supraESR), we introduce a completely new, cost-effective diagnostic parameter, based on an established and easily automated measurement method, that enables low-cost screening for neuroacanthocytosis syndrome, a group of rare neurodegenerative diseases previously detectable only through complex diagnostic tests.

Purpose: Identifying head impacts linked to brain injury in sport remains challenging. Instrumented mouthguards quantify head-impact kinematics, and finite element (FE) modelling can transform these data into brain strain estimates, which may better reflect injury risk than kinematics alone. Here, we examined associations between mouthguard-measured kinematics, FE-derived strain, and plasma brain injury biomarker GFAP following head impacts. Methods: We analysed 41 video-verified impacts from male Australian football players, including 22 assessed for concussion (17 diagnosed) and 19 unassessed. Instrumented mouthguards recorded peak linear acceleration (PLA), peak rotational acceleration, and peak rotational velocity (PRV). Brain strain was estimated using the Imperial College FE brain model, and plasma GFAP was quantified using Simoa. Biomechanical-GFAP associations were examined using Spearman correlations and segmented regression. Results: For impacts overall, plasma GFAP was moderately correlated with PLA ({rho}=0.46, 95% CI: 0.20-0.66), PRV ({rho}=0.53, 95% CI: 0.20-0.78), and strain ({rho}=0.60, 95% CI: 0.32-0.80). Associations were stronger within concussion cases for strain ({rho}=0.86, 95% CI: 0.58-0.97) and PRV ({rho}=0.64, 95% CI: 0.15-0.93). Piecewise regression identified strain levels above which strain-GFAP relationships steepened across the whole-brain and brainstem. In concussion cases, supra-threshold brainstem strain was associated with greater symptoms. Conclusion: Finite element brain strain may better predict brain injury risk following a sport-related head impact than peak acceleration metrics. Stronger associations with plasma GFAP, particularly among concussion cases, and evidence of a biomechanical threshold, support the use of biomarker-informed strain measures in future risk modelling and the development of brain injury screening thresholds.

Mosquito-borne diseases pose a growing public health challenge as climate change reshapes vector population dynamics. West Nile virus (WNV), transmitted between birds and Culex mosquitoes, disproportionately affects Maricopa County, Arizona, one of the nation's highest-burden counties, yet whether models that include weather and avian dynamics improve forecast accuracy remains unclear. Using a 15-year weekly time series of mosquito abundance, mosquito infection prevalence, and human cases, we developed four mechanistic model configurations of varying complexity, from mosquito-human dynamics alone to full models incorporating avian dynamics and weather forcing. We fitted each model to the weekly-observed data, generated probabilistic 1- and 2-week-ahead forecast horizons, and evaluated forecasts against a historical baseline. All configurations fit the data equally regardless of weather or avian dynamics. However, models incorporating both birds and weather created more accurate forecasts of mosquito abundance and mosquito infection prevalence, and all configurations outperformed the baseline for forecasting human cases. Forecast accuracy was highest in summer and fall, and ensemble aggregation sometimes outperformed every individual model, stabilizing predictions across the 15-year record. These findings indicate that avian and weather dynamics are most critical for predicting mosquito-specific data, positioning this framework as a scalable tool for public health planning for WNV surveillance under climate change.

Background National childhood respiratory pathogen spectra are diversifying nearly everywhere - within-country diversity rose in 203 of 204 countries between 1990 and 2023 - yet whether countries are diversifying toward a common spectrum or along divergent paths is unknown. We quantified between-country compositional distance of national pathogen spectra over the same period. Methods We built national pathogen share vectors from Global Burden of Disease Study 2023 lower respiratory infection etiologic attributions (26 pathogens, 204 countries, ages 0-19 years) at five timepoints spanning 1990-2023. Between-country distance was measured as all pairwise Jensen-Shannon divergences (JSD; primary) and Bray-Curtis dissimilarities, with Baselga and Jaccard decompositions; robustness was assessed across metrics, pathogen panels, low-count thresholds and a balanced panel of 107 countries. Results Mean pairwise JSD rose from 0.0084 in 1990 to 0.0283 in 2023 (+238%; trend p = 0.030), peaking in 2021 (+283%) with a partial 2023 pullback. Bray-Curtis dissimilarity rose +120% and the balanced panel +423%. Divergence was entirely balanced variation (share reallocation), with spectrum richness rising from 18.5 to 21.1 of 26 pathogens. Dispersion rose fastest for influenza (coefficient of variation 0.03 to 0.55) and respiratory syncytial virus (0.08 to 0.48). Within-region distance rose in every computable GBD super-region (five of seven): divergence occurs within regions, not between blocs. Conclusions National spectra are re-sorting along country-specific axes as vaccine-preventable dominance recedes at different speeds. Diversification is universal, but convergence is absent: the transition at the etiologic-spectrum level is asynchronous and path-dependent, with implications for empirical treatment policy and pathogen surveillance.

Children and adolescents with hearing loss frequently encounter reduced auditory access and delayed language development, factors that may influence the maturation of executive functions. This study examined developmental differences in planning, a core executive function, in 98 children and adolescents with hearing loss or normal hearing aged 7 to18 years using the Tower of London task. Compared to normal hearing peers, participants with hearing loss made more unnecessary moves and rule violations and initiated problem-solving more rapidly, suggesting reduced preplanning efficiency and increased impulsivity. These group differences were most pronounced in adolescents, who showed faster initiation and greater movement inefficiency than age-matched normal hearing participants. Within the hearing loss group, adolescents displayed higher accuracy and longer initiation times than children, reflecting developmental improvements despite persistent gaps relative to hearing peers. Language development age did not alter the main effects. Findings indicate that reduced early auditory and language access may contribute to differences in planning development, highlighting the need for targeted executive functions support in educational and clinical settings for youth with hearing loss.

Autism Spectrum Disorder (ASD) is a heterogeneous neurodevelopmental condition defined by differences in social communication and restricted, repetitive behaviours. As diagnostic criteria have broadened, ASD is now recognised across a wider range of individuals, raising key questions about its structure: does ASD have discrete sub-types, or is it better conceptualised as a continuous, possibly multidimensional, condition? We aim to explore whether a multidimensional continuum model more accurately captures the variability within ASD. We analysed a large SPARK phenotypic dataset of medical history and diagnostic surveys (background history, SCQ, RBS-R; n=36,710 individuals). We apply and compare two traditional statistical approaches, Factor Analysis and Gaussian Mixture Models, with a modern machine learning technique, the Variational Autoencoder (VAE). VAEs reconstructed unseen test data with ~4-fold better accuracy than Factor Analysis, and ~8-fold better accuracy than Gaussian Mixture Models. We identified four stable latent factors across 100 independently trained VAEs. These four dimensions provide an individual behavioural profile that can be visualized using radar-plots, offering a compact way to compare profiles at the person level. Through further analysis, we found evidence for 3 overlapping clusters or subtypes of ASD identified within the 4D latent space. This work aims to inform new ways of modelling ASD using a VAE that will be able to discern between a continuum or a clustered output and that go beyond binary diagnosis, instead reflecting the complex range of trait profiles, with implications for personalised diagnosis and intervention.

Background: Sickle cell anaemia (SCA) is characterised by chronic haemolysis, inflammation, platelet activation, and recurrent vaso-occlusive complications. Mean platelet volume (MPV) is a readily available platelet index, but evidence regarding its relationship with disease severity in paediatric SCA remains limited and inconsistent, particularly in African populations. Objective: To evaluate the relationship between MPV and disease severity among children with SCA in Kwara State, North-Central Nigeria. Methods: This hospital-based cross-sectional study included 51 clinically stable children with confirmed SCA consecutively recruited from the paediatric haematology clinic of Children Emergency Specialist Hospital, Ilorin. Complete blood count, including MPV, was performed using a Rayto RT-7600 automated haematology analyser. Disease severity was assessed using a composite clinical and laboratory scoring system based on a previously described method. Pearson's correlation, Spearman's rank correlation, simple linear regression, and the Kruskal-Wallis test were used as appropriate. Statistical significance was set at p < 0.05. Results: Of 51 participants, 14 (27.5%) had mild, 33 (64.7%) moderate, and 4 (7.8%) severe disease. Mean MPV was 9.34 +/- 0.76 fL (range, 8.0-11.2). Pearson's correlation showed a weak positive, non-significant linear relationship with severity score (r = 0.231, p = 0.103), whereas Spearman's analysis showed a weak positive monotonic association (rho = 0.286, p = 0.042). Regression explained 5.3% of severity-score variation (R2 = 0.053, p = 0.103). MPV did not differ significantly across severity categories (H = 2.163, p = 0.339). MPV correlated inversely with haemoglobin (r = -0.556, p < 0.001) and positively with platelet count (r = 0.307, p = 0.029). Conclusion: MPV showed a weak relationship with disease severity but inconsistent statistical evidence across analyses. The limited explained variance and absence of significant differences between severity categories do not support MPV as a standalone severity marker. Larger longitudinal studies are warranted. Keywords: Sickle cell anaemia; Mean platelet volume; Disease severity; Platelet indices; Paediatric haematology; Cross-sectional study; Nigeria.

Contact inhibition of proliferation (CIP) enables epithelial tissues to self-regulate growth and maintain tissue homeostasis. However, how cell-level mechanical contact, tissue-scale structural order, and proliferation kinetics interplay remains a fundamental open question in living matter physics. Here, we present a particle-based model of a confluent epithelial monolayer governed by overdamped dynamics, where individual cells interact via a two-dimensional hard core- soft shoulder potential. By comparing structural evolution during quasistatic densification with previously reported experimental division kinetics, we find that the dynamics of proliferation arrest mimics the onset of direct steric contacts between the hard cores of the shell. Identifying hard core contacts as the physical driver of CIP, we couple our mechanical model with a stochastic Monte Carlo division scheme in which the instantaneous division rate decreases to zero from an intrinsic value as the number of hard core contact increases to six from zero. We demonstrate that for high intrinsic division rates, the cellular densification outpaces mechanical relaxation. This kinetic mismatch drives premature hard-core contact formation, shifts the onset of jamming and contact inhibition to lower packing fractions, and induces increasingly disordered transient configurations before the tissue universally converges to a hexagonal close-packed limit. Our model's predicted division kinetics and structural order evolution are consistent with epithelial monolayer experiments, both reported and our own. This minimal physical framework links single-cell steric contact mechanics directly to tissue-scale growth regulation and structural evolution.

Background Sleep, mood, and affective states are mutually connected. There is a paucity of studies, however, that have considered bidirectional relationships between daily sleep-affective dyads in naturalistic settings, particularly for shift workers. Objective To evaluate the dynamic and temporal interplay of daily smartphone-based self-reported sleep measurements, dimensions of affective experience and cognitive processing in UK shift working nurses. Methods The EClocker Study prospectively monitored 102 National Health Service (NHS) nurses (aged 25-61 years, 83.3% female) working standard (day shift) and non-standard (fast rotating shifts) schedules over a two-week period. Smartphone-based Experience Sampling Methodology (ESM) recorded daily sleep, mood, momentary affect and cognitive attentional functioning. Self-reported burnout, emotional dysregulation, emotion reactivity and affective dimensions (positive and negative) were also collected. Findings Overall, NHS nurses reported a high prevalence of depressive symptoms, stress, burnout and sleep-circadian rhythm disturbances. Generalised Additive Modelling (GAMs) revealed that NHS nurses higher perceived sleep quality predicted better next-day mood state, while better daytime mood was associated with reduced sleep onset latency, such that participants reported falling asleep faster. In contrast, daytime mood or affect (positive and negative) had no substantial, direct impact on nurses subjective sleep parameters (sleep quality, sleep duration, sleep efficiency). Exposure to fast rotating night shifts across the two-week study was associated with more frequent response errors on a Choice Reaction Time (CRT) cognitive task, while daytime somnolence did not adversely influence nurses momentary reaction time speeds or attentional function. Conclusions Clinically relevant sleep impairments, insomnia-related symptoms, elevated stress, and poor mood were pervasive in a sample of UK NHS nurses, regardless of shift type. Sleep quality impacted next-day mood and daytime mood impacted sleep latency, while rotating shifts led to an increase in cognitive errors. Recognising the impact of shiftwork and designing interventions to promote better sleep quality offer potential to enhance mood and performance in healthcare professionals. Clinical implications We need to implement and evaluate interventions that regularise sleep patterns and promote sleep quality to alleviate mood symptoms among frontline NHS shift workers.

Digital phenotyping, which is defined as quantifying someone's behaviour with digital devices, provides unprecedented opportunities for understanding human mental health but is hampered by high levels of inter-individual variability. Here, we propose a new method to address this, parsing inter-individual variability by decomposing the digital phenotype dynamics into latent trajectories and using each individual's trajectory membership as a moderator when modelling psychopathology over the same timeframe. We applied our method in the context of mood symptom exacerbation across the menstrual cycle, where symptom severity and timing are inconsistent between individuals. Using the BiAffect platform to collect smartphone typing dynamics, we found stable trajectories in smartphone movement rate: one group of participants showed substantial movement rate fluctuations across the menstrual cycle, whilst the others did not. Participants with movement fluctuations displayed increased fluctuations across the cycle in prospective anhedonia and depression ratings, but not in anxiety, irritability, and suicidal ideation.

Background: Pressure ulcers often develop complex three-dimensional morphologies that extend beyond the visible wound surface. Subsurface extensions such as tunneling and undermining create hidden cavities that complicate clinical assessment and wound management. Despite their clinical relevance, the prevalence and spatial characteristics of these subsurface wound morphologies have not been well characterized at scale. Methods: We performed a registry-based analysis using data from the LIFT-OFF Pressure Ulcer Registry, which captures longitudinal clinical documentation of pressure ulcers treated in routine care. The registry included approximately 18,000 patients with 32,000 documented pressure ulcers. Spatial characteristics of tunneling and undermining were analyzed using measurements recorded during routine wound assessments, including tract length, direction, and circumferential extent. Directional and circumferential distributions of subsurface defects were examined to characterize wound geometry. Results: Tunneling was present in 764 of 14,700 full-thickness pressure ulcers (5.2%), whereas undermining occurred in 2,293 wounds (15.6%). Tunneling tracts were typically short and exhibited directional clustering relative to the wound bed. In contrast, undermining demonstrated broader circumferential distributions and frequently involved larger subsurface separations beneath the wound margin. Both morphologies demonstrated distinct spatial patterns across anatomical locations and wound stages. Conclusion: Tunneling and undermining are common subsurface features of pressure ulcers and exhibit distinct spatial geometries. Whereas tunneling manifests as directional tract-like extensions, undermining more frequently produces circumferential tissue separation beneath wound margins. Improved characterization of subsurface wound architecture may enhance assessment of wound complexity and provide information not captured by surface measurements alone. Future studies should evaluate whether these features contribute to wound severity assessment, prognosis, and risk stratification.

Objectives To analyse the long-term effects of a cruciate ligament (CL) injury on health and socioeconomic outcomes. Methods We used a comprehensive national injury insurance database to identify CL injuries occurring in New Zealand between 2009 and 2022, and employed a doubly robust staggered difference-in-differences research design to identify the effects of these injuries on outcomes up to 10 years after injury. The outcomes of interest were healthcare use (hospitalisations, emergency department visits, medications, knee replacement surgery for osteoarthritis), associated healthcare costs, and labour market outcomes (employment rates, income, and government benefit payments). Results We identified 61 344 CL injuries for inclusion in the analysis. Over 10-year follow-up, a CL injury resulted in increased healthcare use (0.6 more hospitalizations [95%CI 0.4 to 0.7], 1.7 more days spent in hospital [95%CI 1.3 to 2.1], 0.4 more emergency department visits [95%CI 0.3 to 0.6], 2.5 more outpatient visits [95%CI 1.8 to 3.2], and 4.7 more medications dispensed [95%CI -1.8 to 11.2]) and public healthcare costs ($7 537; 95%CI 5 888 to 9 186), reduced income (-$6 060; 95%CI -11 644 to -475), and increased benefit payments ($1 152; 95%CI 542 to 1 761). Conclusion CL injuries have long-term impacts on healthcare use and socioeconomic outcomes. Strategies to reduce the incidence of CL injuries have the potential to realise large health and economic benefits.

Background: Respiratory syncytial virus (RSV) remains a leading cause of hospitalization in infancy, with severe outcomes influenced by both contact patterns and passive immunity. Non-pharmaceutical interventions (NPIs) during the COVID-19 pandemic suppressed RSV circulation and reduced opportunities for maternal immune boosting, potentially altering protection among newborns. We evaluated whether incorporating time-varying maternal immunity improves the ability of an age-structured transmission model to predict post-pandemic RSV hospitalization patterns in infants. Methods: We analyzed population-based RSV hospitalizations among Ontario (Canada) infants (<1 year) from July 2, 2017 to June 25, 2024, using linked administrative databases. A deterministic compartmental model across seven age classes was calibrated against pre-pandemic data using Latin Hypercube Sampling. We compared a model incorporating time-varying contact rates alone against a specification that additionally included time-varying maternal immunity. Results: Both specifications accurately reproduced pre-pandemic seasonality and macro-level post-pandemic resurgence features. The constant maternal immunity model showed slightly better accuracy in capturing the 2021/22 peak compared to the time-varying maternal immunity specification. However, both qualitatively captured the continued near-absence of RSV and the observed peak was captured within the 95% credible intervals. While both models precisely captured the timing and overwhelming surge of admissions that occurred in 2022/23, they failed to capture the premature peak timing and magnitude in 2023/24. Conclusions: Incorporating time-varying maternal immunity did not improve model accuracy post-pandemic. While maternal protection is essential for evaluating infant immunizations, population-level contact shifts primarily shaped post-pandemic RSV seasonality, indicating that models must account for these mechanisms of RSV transmission dynamics.

Generalization is a key cognitive process that allows humans to flexibly apply prior knowledge to guide new behaviors. Difficulties with generalization and flexibility are observed across neurodevelopmental disorders, especially autism, limiting adaptive function and quality of life. Cognitive-behavioral treatment benefits some but not all autistic individuals. As treatment requires application of learned skills to everyday life, variability in generalization ability may limit intervention success in autism. While cognitive substrates of learning and generalization are well established, their potential for explaining clinical outcomes is not known. Here, we combined a category learning task with computational modelling to distinguish two learning strategies underlying generalization -- prototype abstraction vs. exemplar memorization -- and tested whether individual differences in these learning strategies predicted real-world intervention outcomes in autistic youth. Fifty-four participants completed the category learning task at two pre-intervention timepoints, and then completed Unstuck and On Target:14-22 intervention targeting flexible problem solving, goal setting, and planning. We found that participants who consistently relied on prototype abstraction (N=26) were subsequently more likely to benefit from the intervention, showing improvement in parent- and self-reported flexibility. These findings identify prototype abstraction as a clinically relevant cognitive capacity that may help explain individual differences in intervention response and support the tailoring of interventions. More broadly, they demonstrate the value of linking basic cognitive mechanisms to clinical outcomes and may inform strategies to enhance the effectiveness of cognitive-behavioral interventions for youth with developmental disabilities.

Large Language Models (LLMs) are increasingly used by clinicians and patients for medical queries, yet their accuracy and safety at the specialist level in hematology remain insufficiently characterised. We benchmarked ten frontier proprietary and open-weight LLMs across two generations on 1,477 board-style hematology multiple-choice questions (MCQs) derived from five educational datasets spanning nine disease areas and six clinical skill domains, including text-only and multimodal case vignettes. Claude Opus 5 had the highest mean accuracy (92.7% text, 76.9% multimodal), followed closely by Gemini-3.1 Pro (91.4% and 78.7%), Gemini-3.6 Flash (91.0% and 74.8%) and GPT-5.6 Sol (89.9% and 76.7%). Accuracy significantly correlated with model size both for text-only and multimodal MCQs. Between model generations, the largest improvements in accuracy were seen for open-weight models whereas proprietary models showed only marginal gains. In error analysis, top-performing models exhibited highly concordant failure patterns, suggesting shared limitations on challenging cases. Frontier LLMs exhibit substantial specialist hematology knowledge across diverse subspecialist domains and clinical skill sets. Yet, despite high accuracy on board-style questions in hematology, continuous expert-on-the-loop output monitoring is paramount.

A dental AI system and a dentist reading the same radiographs form a paired comparison. Published comparative studies often report the two arms separately against a reference standard, leaving the joint pattern of correctness between them unavailable for secondary paired inference. We show what that omission costs. The accuracy difference remains exactly identified; its sampling variance does not, so the report contains the estimate and not its uncertainty. On a study of 282 units, two published accuracies are consistent with 38 distinct joint tables whose confidence intervals differ in width by a factor of 2.5. The consequence is a three-zone decision map rather than a single threshold: differences at or below 1.06 points are non-significant under every compatible table, differences at or above 6.03 points are significant under every compatible table, and in between the published numbers cannot decide. We then show the omission is repairable at negligible cost. One additional integer, the number of units both arms classify correctly, identifies the joint table exactly and restores standard paired inference. For a panel of readers the pairwise dependences must arise from one joint distribution, a constraint that binds once three readers are present; publishing each reader's joint-correct count against a single reference reader cannot widen and may tighten every pairwise bound, and in a 7-arm experiment reduced them by a median of 37% even for pairs excluding that reference. Where the integer was never published we give DentalPair-Cert, an interval with finite-sample coverage uniformly over every admissible within-unit AI-dentist dependence under the independent-sampling-unit model, certified in both the nuisance maximization and the inversion. Across 4,200,000 simulated comparisons an independence analysis falls to 74.5% coverage with 12.2% type-I error; in a purposive sample of 9 recent comparative studies, 1 reported a paired test on discordant units.

Despite suppressive antiretroviral therapy, many people with HIV (PWH) retain chronic interferon-associated immune dysregulation. Observational data from the Swiss HIV Cohort Study linked asymptomatic mycobacterial exposure to lower viral set points, reduced interferon-associated activity, and attenuated HIV-specific antibody responses, a pattern sharing features with HIV elite controllers and natural hosts of primate lentiviruses. We therefore examined whether Bacillus Calmette-Guerin (BCG) vaccination could induce a related immune configuration in ART-treated PWH. Using longitudinal systems-level profiling within the BELIEVE trial, we found that BCG reduced constitutive NK cell IFN-{gamma} production and PBMC-mediated direct cytotoxicity without impairing inducible cytokine responses or antibody-dependent cellular cytotoxicity. Multiomic and proteomic analyses showed reduced interferon- and activation-associated programs, while adaptive immune parameters remained largely stable and follow-up revealed no obvious adverse clinical pattern. This configuration, reduced baseline interferon activity coexisting with preserved Fc-dependent effector function, shares selected features with immune states described in natural lentiviral control and provides a rationale for testing BCG in combination with antibody-based HIV interventions.

Background Little evidence is available on the epidemiology of different knee injuries at a whole-of-population level. The objective of this article is to provide accurate estimates of knee injury incidence by harnessing the unique comprehensive, population-wide data of New Zealand's universal no-fault injury insurance provider, the Accident Compensation Corporation (ACC). Methods We obtained insurance claims data from ACC covering all knee injury insurance claims approved between 2015 and 2024. We calculated the number of injuries and the incidence rate per 100 000 population, by injury type, year, sex, ethnicity, and age. Results The total number of injuries increased from 184 710 (4 067 per 100 000 population) in 2015 to 244 155 (4 701 per 100 000) in 2024. The most common injuries were other/unspecified ligament sprains, contusions, and collateral ligament sprains. Ligament and cartilage injuries were more common for males than for females, while contusions were more common for females. Ligament tears and dislocations were more common in younger people (15 to 35 years of age), while cartilage injuries were more common at older ages (40 to 65 years). Discussion and Conclusions The rate of knee injuries observed in this study was higher than previously reported in other settings, probably due to broader coverage of injuries treated in primary and community care settings. A broad range of injuries were common, including those that have received less attention in the epidemiological literature to date. More research is needed on the prevention, burden, and outcomes of different knee injuries, beyond a narrow focus on cruciate ligament injuries.

Background: Snakebite envenoming is a major cause of preventable death and disability in the Brazilian Amazon, where long distances, sparse roads, and dependence on river transport delay access to antivenom. We developed location-allocation models to identify community health centers that could strategically expand access to antivenom in Amazonas State, Brazil. Methodology/Principal Findings: We conducted an ecological geospatial study using a 2025 WorldPop population surface, locations of existing and candidate health facilities, and a multimodal road-and-river transportation network derived from OpenStreetMap and HydroSHEDS. Population demand was represented by 7,065 populated centroids, including 1,586 within Indigenous territories. We applied a maximize-coverage algorithm with a six-hour travel-time threshold. Two models were developed: one for Amazonas excluding Manaus and one for populations living in Indigenous territories. Both models began with 77 facilities already providing antivenom and progressively added candidate community health centers until coverage gains plateaued. The plateau occurred at 110 facilities, corresponding to 33 additional centers. In the model excluding Manaus, this configuration covered 1,118,831 people, or 75.11% of the target population; 87.61% of those covered could reach care within three hours. In Indigenous territories, coverage increased from 50.55% to 69.50%, reaching 50,434 people, of whom 81.39% were within three hours of care. Validation used 3,595 snakebite notifications from the 30 highest-burden municipalities in the Brazilian Notifiable Diseases Information System during 2023-2025. The median proportion reaching care within six hours was 40.81% in observed data and 72.17% in model estimates. Conclusions/Significance: Strategically equipping 33 additional existing community health centers could substantially expand timely access to antivenom, particularly in rural and Indigenous areas. Location-allocation modeling that incorporates river transportation can support evidence-based decentralization of time-sensitive health services in geographically complex settings.